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Abstract 

Recent theoretical predictions on DNA mechanical separation induced by 
pulling forces are numerically tested within a model in which self-avoidance 
for DNA strands is fully taken into account. DNA strands are described by 
interacting pairs of self avoiding walks (SAW) which are pulled apart by a 
force applied at the two extremities. The whole phase diagram is spanned 
by extensive Monte Carlo (MC) simulations and the existence of a low tem- 
perature denaturation is confirmed. The basic features of the phase diagram 
and the re-entrant phase boundary are also obtained with a simple heuristic 
argument based on an energy-entropy estimate. 
05.10.Ln, 64.60.Cn, 87.15.La, 87.14.Gy 
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In the last few years, enormous technical improvements of experimental devices such as 
atomic force microscopesS'l, optical tweezersl'i and glass microneedles!!'! have allowed mi- 
cromanipulation of single biological systems, and the determination and characterization of 
new, in part unexpected, behaviours of biomoleculesB. These developments have been par- 
ticularly interesting and innovative for DNA molecules. DNA micromanipulation techniques 
involve handles attached to the two ends of the molecule which serve as anchoring points 
to manipulable physical supports!. In this way it has been possible to study the response 
of DNA to external torquesi0 and its mechanical unzipping in the absence of enzymes0'0. 
Nevertheless, only very recently theoretical investigations have started to consider the denat- 
uration of DNA under the presence of forces or torquesEHll. On the other hand theoretical 
studies of thermal denaturation (or melting) of DNA have a much longer history§'§(for 
more recent studies seeEUli) 

Recently, Bhattacharjee@ has extended a minimal model@ in which the two strands of 
DNA are ideal polymer chains, introducing a pulling force applied to two extremities^'^. By 
mapping into a non-hermitian quantum mechanics problem he has shown that the relative 
polymer can be pulled out (i.e. the DNA can be unzipped) only if the applied force exceeds 
a critical value. 

These results have been further developed^S by considering analytically solvable models 
in which DNA strands are represented by directed walks. Surprisingly the force vs temper- 
ature phase diagram revealed the presence§!'@ of a novel re-entrant unzipping transition at 
low temperature (see Fig. 1). These results have been subsequently confirmed by other 
authonS@ within ideal chain models of DNA. 

The aim of our paper is to study the phase diagram of DNA in the presence of a pulling 
force with a more realistic model in which self-avoidance for the DNA strands is explicitly 
taken into account to verify the robustness of the low temperature re-entrance. 

We modeJ^E the two strands of a 2N base-pair DNA by two interacting iV steps SAW, 
on a cubic lattice, T 1 = {Tq, T\, . . . , Tjy} and T 2 = {Tq, 1^, . . . , T 2 N }, with T l p as the lattice 
site occupied by the pth monomer of the zth chain (i = 1,2). The first monomers T^Tq 



are kept fixed at one lattice unit apart. Each monomer (site) on a strand corresponds to a 
base and it is supposed to have its complementary at the same contour position in the other 
strand. Two complementary monomers are favoured to be in contact (i. e., if they occupy 
two nearest neighbour sites of the lattice) by an energetic gain e representing the binding 
energy. 

In this approach0El~ii] the sequence of bases is not explicitly considered since the model 
is coarse grained in character. In this respect each monomer is not to be thought of as a 
single base, but as a group of bases (block). Consequently a mismatch between correspond- 
ing blocks has to be very disfavoured with respect to a correct matching. This aspect is 
adequately treated by our model. 

To describe the mechanical unzipping of DNA an energetic term — / • r is added where 
r = rjy — T 2 N is the vector joining the two ends of the strands and / is an external force 
taken, without loss of generality, along the z axis, e.g. / = fe z (see Fig. 2). Note that at 
/ = 0, this model is in the same universality class as that of Ref.0. 

The Hamiltonian of the problem is then given by: 



where 8 (a, b) is the Kronecker delta and \x — y\ is the Euclidean distance between x and y 
in lattice units. 

Throughout all the calculations e = +1 and therefore all the thermodynamic properties of 
the systems depend on the temperature T = and on the force /. The free energy per 
monomer .Fjv(/3, /), is 



where the sum is over all the possible pairs of N-step SAW's T 1 and T 2 . 

The phase diagram can be studied by looking at the average fraction number of contact 
and at the z-component of the average end point separation (r z ) defined as: 




(1) 




(2) 



In the denaturated regime should vanish whereas it should be nonzero in the zipped 
phase. To compute averages, for fixed values of / and N, Monte-Carlo simulations of the 
model have been performed by a Metropolis based hybrid algorithm0S that consists of local 
and global (pivot moves) deformations attempted on each strand of the polymer. In addition 
a move that tries to interchange the position of pairs of zipped (or unzipped) portions of 
the double chain has been considered. 

The efficiency of the algorithm has been further enhanced by the multiple Markov chain 
sampling schemeHS, a method that recently has been shown to be quite effective in ex- 
ploring the low temperature phase diagrams of interacting polymers. All the runs have 
been performed by covering, with at most 30 different Markov chains, a large portion of the 
temperature space ranging from T = 3 down to T = 0.1. The values of N considered range 
from 2N = 50 up to 2N = 150 and the values of / go from / = up to / = 1. 

In Fig. 2 and Fig. 3 the numerical results respectively for and (r z )/(2N) are plotted 
as a function of T for different values of the force strength / and for systems of 2N = 150 
monomers. The data signal clearly the presence of three distinct regimes. 

a) For / < /o ~ 0.5 the system undergoes a transition from the denaturated to the 
zipped regime. In particular at / = results similar to those found in ref.0 are re-obtained, 
namely a first order like transition at To = 0.58 ± 0.02, where however the fluctuation in 
the size of the "bubbles" diverges. Interestingly, our data support the prediction that the 
critical force near the melting point (Xb) is / ~ \To — T\ a with a < lHH (see also below). 
This is different from the 3-dimensional Gaussian chain caseS or the directed case of Fig. 
lH (in d — 1 + 1) where one finds a = 1. In these two cases the zero- force melting transition 
at To is continuous. 

b) For /o < / < ft ~ 0.70 the unzipped/zipped transition is still present but if the 
temperature is further lowered the two strands separate again through a low temperature 
denaturation or unzipping induced by the presence of the force. This is established by the 
peaks in Fig. 3 and by the minima in Fig. 4 for / = 0.625 and / = 0.65. 

Analysis of the specific heat peaks for different chain lengths support the fact that both 
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these transitions are first order. One might think of the existence of a subtle difference 
between the low and the high temperature denaturation. To explore this possibility we 
studied the size of the molten "bubbles". For / = a simple argument^ predicts that 
the size fluctuation of this "bubble" diverges when the thermal denaturation temperature 
Tjy is approached from below. As already mentioned we observed this effect for / = 0, 
but not for / 7^ 0, both approaching the low temperature transition from above or the 
high temperature transition from below. We explain this fact by noticing that for / 7^ 
the transition temperature is less than To and the part of the two strands containing the 
"bubbles" is not influenced by the presence of the external force (the external force act as a 
boundary effect). Thus, by approaching the transition line from the zipped phase, bubbles 
can never become critical if / 7^ 0. Therefore we conclude that from a statistical point of 
view low and high temperature denaturation transition are equivalent. 

c) For f > ft the pulling force is always dominating over the base-base interaction and 
the system is always in the denaturated state. 

The value of /o can be exactly calculated and the basic features of the phase diagram 
and the re-entrant phase boundary can be obtained with a simple heuristic argument based 
on an energy-entropy estimate. Indeed at T = 0, when thermal fluctuations are absent, 
the only allowed configurations are those in which the first I monomers of the strands are 
parallel and form I contacts, whereas the remaining N — I are pointing in opposite directions 
parallel to the force. Hence, the resulting energy is a simple function of I: 

E{l) = -l(l-2f)-f(2N + l) (4) 

For / > 1/2 the energy minimum occurs for I = (complete stretching) whereas for / < 1/2 
it occurs for I = N (complete zipping). This argument states definitively that /o = 1/2 in full 
agreement with our numerical findings. Notice that for / = f = 1/2, E(l) is independent 
of I, indicating coexistence. 

For low enough temperatures we expect that the dominant configurations are still of 
F-like configurations with 2(iV — I) monomers zipped and two "free" strands each of length 
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I whose extrema are kept at a relative position r by the force /. In this Y model the energy 
is thus given by: 



E = -(N-l)e-f- f 



(5) 



whereas the entropy is 



S = lot 



21..N-1 



(6) 



where P/v(r) is the probability that an iV-step strand has an end-to-end displacement r 
and [i u and \x z are the effective coordination of the single strand and of the double strand 
respectively. For the simple model introduced before \i z = fi u but in general this need not 
be true. 

It is well knownS that 

s 



Pn{r) ~ exp 



r «r«l 



(7) 



where r = \r\, 5 = (1 — v is the usual polymer critical exponential and c is a constant. 
For an ideal polymer v — 1/2 and one recovers the usual Gaussian distribution. 

In the last equation power law corrections lead to subleading terms in the entropy esti- 
mate. From the last three equations the following free energy is obtained: 

s 



2 

F(l,r) =l{e-T\og^L) -r ■ f + cT (^) + <-onst 

fJi Z \L 



(8) 



Minimization of F with respect to r yields 



r(I, /) = I 



c5T 



(9) 



where r is the magnitude of r in the direction of /. This r gives a free energy 

A 8-1 



F(l)=F(l,r(l,f)) = l 



e — T los 



5 s - 1 



(10) 



implying that F(l) gets its minimum at I = (zipped phase) and / = N (unzipped phase) 
when / < f c (T) and / > f c (T) respectively with 
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/cCO~t 1 -"(i-^) i ' ) (11) 

To = e/logf^"-) being the thermal denaturation temperature for this Y model. Thus at 
f c (T) a first order phase transition occurs where the fraction of the bounded bases has a 
jump. The last equation predicts a re-entrance with f c — > when T — > 0. This limit is given 
by a spurious effect due to the assumption that Eq. |7| applies also at very low T where eq. 
H predicts r — > oo. Indeed if Eq. [7| is assumed to hold till P2z( r )/ i « ~ 1> e -§- r — ~ 2/, 
and after that r remains constant, then we easily find that at low T 

f c (T) = ^(e + T log fi z ). (12) 

Eq. fllTD near the denaturation temperature and eq. (O) at the lowest T predict a behaviour 
qualitatively in agreement with the one described in fig. 10 and confirmed by our simulation 
analysis. 

In conclusion we have shown how a pulling force applied to the extrema of a model of 
double stranded DNA molecule could introduce non trivial effects in the phase diagram. 
The results obtained on MC simulations in which, unlike previous studies, excluded volume 
effects are fully taken into account, confirm the existence of a low temperature denaturation. 
Although in our approach the intrinsic helicity of the DNA and the randomness of the bases 
sequences have not been incorporated, we believe , on the basis of preliminary numerical 
results, that they should not affect the main features of the phase diagram and in particular 
the existence of a re-entrance region. Therefore we expect that such an effect could be 
pinpointed by high precision measurements of DNA denaturation under a pulling force. 
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FIG. 1. The exact force (f) versus temperature (T) phase diagram as obtained in ref.E3 for a 
directed model (d=2) of DNA. See eq. (1) for definition of force. The values of ft, fo and Tjj are 
model dependent. 
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FIG. 2. A typical configuration of a double stranded DNA chains (2iV = 200) obtained with 
the model described in the text. The different colours indicate the two different strands. 
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FIG. 3. Average fraction of bound pairs versus T, for several values of the force / and for 
2N = 150. Error bars are of the size of the symbols. The two peaks (/ = 0.625 and / = 0.65) 
signal a zipped phase between two unzipped regions (re-entrance). 
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FIG. 4. End to end end separation in the force direction versus T, for different values of the 
force / for 2iV = 150. Error bars correspond to one standard deviation. The two minima (/ = 0.625 
and / = 0.65) signal a zipped phase between two unzipped regions (re-entrance). 
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